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ABSTRACT 



A two dimensional model of an anchor cable subjected to 
a uniform current, fixed at the sea bed and moving with a 
specified motion at the surface is developed. 

A steady state solution is developed and used as an 
initial condition for the numerical solution of the dynamic 
equations. The method of characteristics is used to carry 
out the numerical integration. 
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TABLE OF SYMBOLS 



A = cross sectional area (in^) 

A n ' = cable acceleration in the normal direction (ft/sec^) 

A s ' = cable acceleration in the tangential direction (ft/sec2) 
C 3 = drag coefficient for a cylinder 
d = diameter of cable (in) 

D' = drag force per unit length of cable (lbs/ft) 

Yg 0 ' = depth measured from the mean water line (ft) 

E = modulus of elasticity (lbs/in^) 

K = cable spring constant (lbs) 

L = length of cable (ft) 

M = mass per unit length of cable (lb-sec^/ft^) 

M' = virtual mass per unit length of cable (lb-sec^/ft^) 
p = density of seawater (lb-sec^/ft^ ) 

¥ = angle between the horizontal and a tangent to the 

cable at a point on the cable 

S' = length along the cable measured from the fixed point 
on the sea bed (ft) 

T' = tension (lbs) 

tg = time scale (period for one cycle) (sec) 

T Q = steady state tension of the cable at the surface (lbs) 
t' = time (sec) 

V Q = velocity scale (velocity amplitude at the surface) 
(ft/sec) 

V' = current velocity (ft/sec) 

V s ' = tangential cable velocity (ft/sec) 

V n ' = normal cable velocity (ft/sec) 
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W = wieght per unit length of cable (lbs/ft) 

Y s ' = vertical distance between the fixed point and the 
top of the cable (ft) 
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TABLE OF DIMENSIONLESS PARAMETERS 



V n = V n '/V 0 

V s - v s -/v 0 

T = T ' /T 0 
A = gt 0 /v 0 
B = T 0 t 0 /LMV 0 
C = t 0 W/V 0 M' 

D = V 0 pC d dt 0 /2M* 
E = V/V Q 
F = L/V 0 t 0 
R = K/WL 
C s = (RA/F) 3 * 

C n = (BCT/AF) ^ 

Ys 0 =^ So '/L 
Y s = Y s '/L 
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I. 



INTRODUCTION 



The study of anchor cable dynamics in recent years has 
been of some importance to those interested in the motion 
of vessels in seaways. The mooring of buoys for navigation 
or the mooring of large tanks to hold oil or other sea- 
transported goods are further examples of the importance of 
this study. The motion of ships riding at anchor in limited 
areas would be of interest, for example, to shipping compa- 
nies that service small ports. 

In a survey [Ref. 1] done by Casarella and Parsons, it 
was mentioned that the first experimental work on mooring 
cables was probably done in 1917 by Relf and Powell. Since 
then several people have worked on the mathematical model 
of this problem. 

Through the years the model has changed to some extent. 
Some of the extensions considered have been the inclusion 
of tangential drag on the cable, the effect of surface waves, 
the effect of current, the effect of added mass, and also 
the effects of cross currents in the three dimensional 
models [Ref. l]. Most of the studies done have been carried 
out for the static configuration; however, dynamic studies 
.with arbitrary motion prescribed at the surface end of the 
cable have also been carried out. 

With the advent of fast computers, time required to get 
numerical solutions has been greatly reduced; however, 
computer time is still considerable as will be shown in this 
analysis . 
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Wilson [Ref. 2] arrived at a model very similar to the 
one used in this thesis. Wilson considered the effect of 
surface waves as well as tangential drag on the cable; 
however, the computer program used to solve the problem was 
proprietary to the company which supported his research and 
was not published. 

For the present analysis, the drag normal to the anchor 
cable only is considered and the effects of waves are 
neglected. The assumption of zero tangential drag has been 
included in the general approach to the problem and seems 
to be valid for reasonable currents and cable configurations. 
Surface waves are also neglected in this analysis as far as 
drag is concerned. The current applied is steady, uniform 
and in the horizontal plane. 

In the dynamic study, the method of characteristics is 
used to solve the equations of motion using a digital computer. 
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II. THEORY 



The problem under consideration is depicted in Figure 1. 

A cable of diameter d and weight W pounds per foot is anchored 
on the sea floor and is suspended from a point on the surface 
as indicated in Figure 1. A uniform current of magnitude V' 
is directed in the positive X-direction and the motion of 
the upper end of the cable may be specified in some suitable 
manner. 

Writing Newton's Law for an elemental length of cable 
as depicted in the insert of Figure 1, the following dynamic 
equation is obtained: 

T'AY - WAs' cos Y - D'As' = M'A n 'As' (1) 

where A n ' is the normal acceleration. 

Dividing equation (1) by As', and taking the limit as As'-^-O 
yields : 

T' — — - W cos Y - D' = M'A n ' (2) 

9s' 

The term M' represents the virtual mass of the cable per unit 
length and D' is defined as the cable drag per unit of length. 
T' denotes the tension in the cable and the angle Y denotes 
the angle between the tangent at a point on the cable and the 
horizontal. In summing the forces in the tangential direc- 
tion the added mass is taken as zero and as a result the 
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FIGURE 1 DEFINITION SKETCH 
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actual mass is simply the cable mass. Also, the tangential 
drag is omitted as stated in the introduction. From the 
forces shown in the insert of Figure 1, the equation of 
motion for the tangential direction takes the form: 

-T ' + (T ' + ill As ' ) cos A4> - WAs ' sin (t + 11) = MA 'As' (3) 

3 s 2 s 

Again, dividing through by As' and taking the limit as As'-*0 
yields : 

|11 -■ W sin Y = MA s ' (4) 

It is necessary also to develop expressions for the 
normal and tangential accelerations in terms of the normal 
and tangential velocities. Referring to Figure 2, the normal 
and tangential acceleration components as well as the corre- 
sponding velocity components can be written in terms of the 
x and y components as 
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in which X, Y, X, Y, denote the velocities and accelerations 
in the X and y directions, respectively. Taking the partial 
derivatives of equations (7) and (8) with respect to time 
and using the definitions of V ' and V g ' as given in equations 
(7) and (8) it is possible to show that the normal and 
tangential accelerations may be written as: 
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The primes indicate that these terms are dimensional terras. 
Dimensionless forms of these parameters will be derived later. 

Substituting equations (S) and (10) for the acceleration 
into equations (2) and (4) , the following expressions for the 
two equations of motion are obtained: 

u/ S V 1 ^ ur 

T' W cos V - D' = M' ( 31 _ + V ' _) (11) 

3s' 3 1 1 S 3 1 1 

~ - W sin t = M - V n ’ 1L_) (12) 

3s ' 3 1 1 3 1 1 



y 




FIGURE 2 Coordinate Systems 



From the dynamic equations, (11) and (12), it can be seen 
that there are four dependent variables, V 1 , V s ' , T',T , all 
of which are functions of s' and t*. In order to solve the 
equations of motion two other equations are necessary. These 
are the two kinematic relations whose derivations follow. 

With the assumption of a linear relationship between the 
stress and strain and referring to Fig. 3, it is possible to 
write the following relationship between the rate of increase 
of tension and the rate of strain, 

T' + 111 As' 




FIGURE 3 Tangential Forces on an Element of Cable 



= -1- [(V '+ j Vs — As') cos AT - V s ' - (V n '+ As') 

At' As' s 9s' s n 3S , 

sin AT] (13) 



which, upon taking the limit as As'->0, becomes: 

1 _ 3 T 1 _ 3 V s ' i 3 t 

K 3 t 1 3s' n 3s' 



(14) 
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Equation (14) states that the rate of change of tension 
equals the "spring constant" times the rate of stretching. 

In the case of a solid cable of uniform cross section, K is 
simply AE where A is the cross sectional area and E denotes 
the modulus of elasticity. In the case of, for example, a 
chain or stranded cable the value of K would become somewhat 
more complex and would most likely be determined experimentally. 





As 1 



FIGURE 4 Cable Kinematics 



A second kinematic relationship may also be introduced 
in the following manner. The angular velocity of a segment 
As' of the cable may be specified in terms of the velocity 
of its end points. Referring to Figure 4, the angular 



velocity may be written as 

V n 1 + IV 
8T _ n 



AS' - V ' 

8 s 1 n 



+ V< 



, ST As' 



at' 



AS ' 



8S ' As 1 



(15) 



which becomes upon taking the limit as As'->0, 
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The four equations which have been developed, equations 
(11) , (12) , (14) \ (16) , represent the complete system and 

are all that are required to determine the four dependent 
variables V n ' , V g 1 , T * r T , in terms of the functions s' and 
t'. In summary, the basic system of equations is: 



t ' ii 
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The four equations which have now been specified equa- 
tions (17) through (20) , may be integrated numerically to 
obtain a solution. It is possible to carry out this numerical 
integration using a finite difference formulation in the 
s'-t' plane with a rectangular grid. However, it is also 
possible, and more convenient, to use the method of charac- 
teristics as outlined by, for example, Crandall [Ref. 3]. 

Using his suggested method it can be shown that the system 
is hyperbolic and equations (17) through (20) can be trans- 
formed into the form of first order differential equations 
valid along their respective characteristics. 

According to Crandall, four other relationships in 
addition to equations (17) through (20) are needed to carry 
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out the integration by use of the method of characteristics. 

These relations are directional derivatives for the dependent 

variables V_ ' , V c 1 , T' and ¥ as defined by 
n fa 
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Equations (21) through (24) denote incremental changes 
associated with small displacements ds ' and dt ' in the 
t'-s' plane. In order to determine the characteristics of 
the set of equations, equations (17) through (24) are set 
in a matrix form as follows: 
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Expanding the determinant of the square matrix and 
equating it to zero the following result is obtained: 



K , T ' x j , 2 ^ 4 . . 2 A T'K **..4 



ds ' 4 - ( * + ^_)ds l2 dt 
M M' 



+ dt' 

MM' 



= 0 



( 26 ) 



Dividing by dt 1 4 equation (26) becomes: 



(Sli) 4 - < E + II )(Ell) 2 + I^= o 



dt ' 



M M' dt' MM' 



(27) 



Equation (27) is a quartic, the solution to which will 
produce four real characteristic directions. 

Letting 
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and using these definitions equation (27) becomes: 



(^£i ) 4 - (c • 
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(30) 



Solving now for ds'/dlt,' gives the four solutions: 
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These slopes define the four characteristics in the s'-t' 
plane associated with the system. The values C s ' and C n ' 
agree with those of Wilson [Ref. 2], The fact that there 
are four real characteristics classifies the system as 
hyperbolic . 

Having now determined the slopes of the characteristics, 
the next step is to determine the first order ordinary 
differential equations along each of these characteristics. 
To do this, the right hand side of the system matrix 
equation (Eq. 25) is substituted into any column of the 
square matrix and this determinant is expanded and set 
equal to zero. The first determinant is given in eq. (33) . 
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Expanding this determinant the 



following result is obtained: 
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where C 1 , and C ' are defined by equations (28) and (29) 

II b 

and C = ^s ' . 

dt ' 

It may be noted that when C = ±C n ' both sides of the equation 
vanish. When C g ' is substituted for C the following two 
equations, valid along the ±C S 1 characteristic, are obtained: 



(+c s '), 



D(Vs l ) 
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+ V 1 
v n 



D( y ) 1 D(T') 
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(36) 



Equations (35) and (36) are the ordinary differential 
equations valid along the ±C S ' characteristics; one for the 
positive slope and one for the negative slope. It is 
interesting to note here that the C s ' slopes are constant 
and therefore do not vary throughout the complete s'-t' 
plane . 

Equations for the C n ' slopes are determined in a similar 
manner by substituting the right hand side of the matrix 
equation system (Eq. (25)) into another position in the 
square matrix; the fifth column will produce the desired 
result. After expanding, the following equation is obtained: 



(C s ,2 -C 2 ) 



+ C(C 2 -C 2) EfEl + V (C s '2-c2) 

Dt' S Dt' 



D(¥) 
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M ' 



( 37 ) 
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In this case, if C is replaced by ±C S ' the equation 
becomes trivial. However, if ±C n ' is substituted for C 
the following equations, valid along the second character- 
istics C n ', are obtained: 



(+c n ) 



( - C n ) 



D ( V n 1 ) 
Dt' 

D (V n 1 ) 

_ + 

Dt ' 




DT 

Dt 7 ’ 



(V s ' + 



C ' ) 
n ’ 



DT 

Dt’ 



W cos T + p 1 
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^ W cos T + D' 

M' 



(38) 

(39) 



Equations (35) , (36) , (38) , and (39) represent four 

ordinary first order differential equations which are valid 
along their respective characteristics. The total deriva- 
tives 5 — are defined according to 2 1 ±C j* , 

Dt' Dt' 9 1 1 as' 

ds 1 

where C = denotes the slope of the characteristic in 

the s'-t' plane. Since these equations are first order 
equations they are easily integrated along the characteristic 
directions . 

With a uniform current directed as shown in Figure 1, 
the drag per unit length acting normal to the cable may be 
written in terms of the fluid velocity relative to the cable 
in the form 



D' = 33pC d d(V’ sin T + V n ')|V' sin T + V R ' | (40) 

where p denotes the fluid density, denotes the drag 
coefficient, V denotes the current magnitude and d denotes 
the cable diameter or frontal area per unit length. The 
absolute value is used in order to maintain the sign of the 
relative velocity. 
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A. DIMENSIONLESS PARAMETERS 

The four characteristic equations (35) , (36) , (38) , and 

(39) can now be put in a dimensionless form. Various terms 
need to be defined before the actual analysis is done. 

The solution to the problem as depicted in Figure 1 
depends on a number of parameters. That is, the periodic 
motion at the surface may be characterized by a period (t ) 
and a velocity amplitude (V D ) . Also, at steady state, the 
tension at the surface T Q as well as an angle between the 
horizontal and a tangent Y to the cable may be specified. 
The unstretched length (L) of the cable can be used to 
normalize all length scales and the velocity amplitude may 
be used to normalize the velocities. Applying these 
definitions to the characteristic equations the following 
dimensionless parameters may be defined: 
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( 47 ) 
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Applying the dimensionless groupings (41) through (56) 
to equations (35) , (36) , (38) , and (39) , the characteristic 

equations can be written in the following dimensionless form 



(+c s ) - D(Vs> + V £1IL + 2_ 5M 

Dt n Dt C s Dt 

l-r \ D(Vs) D(y) , B_ D (T) 

S Dt " n Dt C s Dt 
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-A sin V 



(57) 

(58) 
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— + (V s -FC n ) — = - C cos V - D[E sin S' + V n ] |E sin V + V | 

Dt Dt 



(59) 




Equations (57) through (60) may now be solved numerically 
on the digital computer to obtain the dynamic solution. 
However , in order to proceed with the solution it is first 
necessary to determine an initial condition as well as 
appropriate boundary conditions on the surface and ocean 
bottom. 

B. STEADY STATE SOLUTION 

The steady state solution may be obtained by setting V n , 

V s and 1 equal to zero in equations (57) through (60) . 

9 1 

Accordingly, the following steady state form of equations 
(57) and (58) is obtained: 



which is equivalent to equation (12) if the parameter defini- 
tions are substituted. Equations (59) and (60) become: 

FC n ^ = C cos Y + DE^sin^ Y (62) 

ds 

Equation (62) is equivalent to the equation of motion for 
the normal direction, equation (11). Since the characteris- 
tic equations agree with the equations of motion as far as 
steady state is concerned, equations (59) and (60) will be 



b 

ds 



A sin V 



(61) 
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the governing equations for the steady state solution. 

Also, the fact that both sets of equations agree is an 
indication that the characteristic equations were derived 
correctly . 

Substituting equation (56) for C n and dividing equation 
(61) by equation (62) provides the following expression: 



dT sin ¥ d¥ 

— ; 2 O 

T cos ¥ + DE sin^ ¥ 

C 



(63) 



This equation can be integrated analytically between the 
limits of T to T Q and ¥ to ¥ . Upon integration the follow- 
ing result is obtained: 




(64) 



If there should be zero current then equation (63) would 
reduce considerably since the parameter E would drop out. 
An analysis is done in Appendix B to compare equation (64) 
without the current term to the solution obtained directly 
from equation (63) setting E equal to zero. 

Once T is known from equation (64) then the following 
expression can be used to get s: 



ds 



BCT d¥ 

AC cos ¥ + ADE 2 sin 2 ¥ 



(65) 
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Equation (65) can not be analytically evaluated, and there- 
fore it was integrated numerically on bhe digital computer. 

C. BOUNDARY CONDITIONS 

With the steady state solution used for the initial 
values only boundary conditions at the surface and sea bed 
need be specified. At the sea bed the normal and tangential 
velocities will be zero since the cable is anchored. With 
these two relations plus equations (58) and (60) , the four 
dependent variables can be determined along the sea bed. 

At the surface, equations (57) and (59) can be used; 
however, two other equations need to be specified such that 
the four dependent variables can be calculated. These two 
relations can be specified in several ways. For example, 
the surface velocities may be written as: 
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ii 
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* + V x 


cos 
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ii 


cos 
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< 


sin 
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(67) 



where V x and Vy may be specified functions of time. 

The two relationships required at. the free surface may 
also take the form of the two-dimensional equations of 
motion for a floating vessel. In this case tension would 
be related to the position, velocity and acceleration of 
the floating object. 
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III. NUMERICAL PROCEDURE 



Equations (57) through (60) represent four ordinary first 
order differential equations valid along their respective 
characteristics. In order to carry out the numerical integra- 
tion, these equations are placed in the following difference 
form: 




( 68 ) 
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Figure 5 System of Characteristics 
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The subscripts 1, 2,3,4 denote the initial point where it 
is supposed that all quantities V g , V , T, and t are known. 

In order to "march out" and determine new values at the new 
location denoted by the subscript 5, the four linear equations 
(68) through (71) are solved simultaneously. Letting the 
right hand sides of equations (68) through (71) be equal to 
AA, CC, DD, and BB, respectively, the resulting expressions 
for the dependent variables at the new location 5 are: 
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The basic method of calculation indicated in equations 
(72) through (75) was used to determine values of the depen- 
dent variables V , V , T, and Y throughout the t-s plane. 

The numerical scheme utilized consisted of first dividing the 

s axis into equal subdivisions. The solution for V , V , T, 

s ri 

and Y at t=0, or initial condition, was obtained from the 
steady state solution of equations (64) and (65) . The steady 
state solution provided values for V , V , T, and Y at all of 
the nodal points along the S axis. The layout in the t-s 
plane is shown in Figure 6 . 

Using equations (72) through (75) a new set of values 
were then determined at the next row of points located at 
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time Dt = Ds/C from the fe axis. Once all the values of the 
s 

four dependent variables were determined at this value of 
time, the solution was marched out another time step, Dt, so 
that values of the dependent variables were eventually 
obtained for a large area on the t-s plane. 

It may be noted in Figures 5 and 6 that the ±C n charac- 
teristics do not originate from a nodal point where values 
of the dependent variables are calculated. In the numerical 
procedure the method of calculation involved a linear inter- 
polation between points 1 and 4 in Figure 5 in order to 
evaluate the dependent variables at points 2 and 3. 

As indicated previously, the steady state solution was 
utilized to determine the initial values of the four dependent 
variables along the s axis. On the t axis, corresponding to 
the ocean floor, it is also necessary to apply some kind of 
boundary condition. In the example worked out, this boundary 

condition was taken as V = V =0 which corresponds to a 

ns 

fixed anchor point. 
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FIGURE 6 Layout of Characteristics 
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IV. DISCUSSION OF RESULTS 



In order to work a numerical example, it is necessary to 

specify values for the parameters appearing in equations (41) 

through (56) . For the example worked, a steel cable of uniform 

cross section with a diameter (d) of 2 inches and a modulus 

of elasticity (E) of 30 x 10^ psi was selected. The density 

of the steel was 480 lbs/ft"^ and the cable was 1000 feet long. 

The period (t ) of the motion at the surface was 10 seconds 

and the velocity amplitude V q was 3.0 feet/second. The 

2 4 

density of seawater was taken to be 2.0 lbs-sec /ft . The 

. ird^ 

virtual mass as described by M'=M + C p — * — where C denotes 

J nr 4 m 

the added mass coefficient [Ref. 2], and p the density of 
seawater, was . calculated to be 3.49. The added mass coeffi- 
cient used in the above calculation was taken as 1.5. The 
drag coefficient was taken as 1.5 which, according to 
Wilson [Ref. 2], is valid for Reynolds numbers less than 10 . 

It is now possible to calculate values for the parameters in 
equations (4 8) through (56) . 

With an accurate steady state solution for tension having 
been calculated by equation (64) , a check on the numerical 
integrations done in the characteristics method is possible. 
That is, it is possible to apply steady state boundary condi- 
tions at the free surface and recompute the steady state 
solution using the method of characteristic procedure. These 
results can then be compared with the original steady state 
solution. 
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In order to do this the steady state boundary condition at 
the surface was written in the following manner: 

T cos y| =cos Y (76) 

S=1 o 

This equation states that the horizontal component of tension 
at the surface was held constant. The vertical component was 
written as: 



T sin 4*1 . = sin Y + C (Y - Y ) (77) 

1 s=l o o s s Q 

Equation (77) is a simplified expression representing the 
vertical component of tension as caused by a massless vessel 
floating on a still water surface in water of depth, Y 

o 

The dimensionless water depth Y is defined as Y = 

. s s 

\ o o 

/' sin Y ds | _ and likewise the vertical span of the cable 

0 ' 1 

at time t is defined as y = / sin Y ds| . The product 

s • 0 

C (y - y ) represents a buoyant force caused by a sinking 
° o 

of the vessel a distance of (y - y ) . The constant C is 

o 

used to describe the interaction of the ship with the free 
surface. The "spring rate" C Q represents the slope of the 
buoyant force-sinkage curve and is accordingly defined as 
C = (A . ) y ) L/T where A . is the water line area, y 
is the specific weight of seawater, and L is the length of 
the cable. Using equations (76) and (77) as the surface 
boundary condition, the method of characteristics was "marched 
out" several time steps in order to compare the results 
obtained from the method of characteristics procedure with the 
initial steady state results. 
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The results of the above described procedure are shown 
in Figure 7. The results corresponding to t=0 were calculated 
by use of equations (64) and (65) . Using these initial values, 
the method of characteristics procedure was then used to 
continue this steady state solution for some time. It may be 
seen from Figure 7 that some drift in the results took place 
due to the numerical integration by the method of characteris- 
tics when compared to the initial solution. The figure shows 
that the vessel sank slightly as the tension increased. This 
change from the initial solution is not to be construed as 
inaccuracy, however. The new steady state solution simply 
represents a second possible steady state solution. 

Since the steady state solution using the method of 
characteristics apparently provided good results, it was 
decided to add a dynamic boundary condition at the surface in 
order to induce a dynamic response of the cable. The boundary 
conditions at the surface, equations (76) and (77) , were 
therefore modified slightly to include motion. Equation (76) 
was applied so that the horizontal component of tension was 
held constant. The vertical component of tension, however, 
was expressed in the following way: 

T sin ? I , = sin v + C [ (Y -Y ) + A sin 2irt] (78) 

' s=l o o'- s s m 

o 

Figure 8 depicts the physical problem defined by the 
boundary condition equation (78) . The equation describes the 
vertical component of tension as caused by a massless vessel 
floating in a wave of elevation A sin 2ut measured with 
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FIGURE 7 Steady State Results 
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FIGURE 8 Dynamic System 
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respect to the mean water level. The term A^ represents the 
wave amplitude and is made dimensionless with the cable length. 
The terms within the square brackets represent the displace- 
ment of the vessel from its mean position with respect to the 
free water surface. As defined previously, C Q denotes the 
slope of the buoyant force versus sinkage curve. 

Figure 9 is a graph of the dynamic response of the cable 

to the above described motion. On this graph are plotted the 

dimensionless tension and depth against a dimensionless time 

scale. Tension is made dimensionless with T which is 

o 

defined as T q = WL/sin ¥ , ¥ = 45° being the initial steady 

state angle at the surface. The depth is made dimensionless 
with the cable length, which is 1000 feet for the example 
worked. Finally, the time is made dimensionless with the 
period of the motion at the surface which is 10 seconds. The 
T q calculated from the above expression is approximately 
14,150 pounds. The results of the steady state analysis. 

Figure 7 is also shown on Figure 9. Figure 9 shows a sinu- 
soidal response which is expected from the boundary condition 
driving the motion. Although the average tension at the 
surface should be near 1.0, it can be noticed that after 
about three cycles the tension appears to be oscillating 
around a mean value which has a 6 per cent shift. This new 
mean position is equal to the one calculated from the method 
of characteristics under steady state conditions. 

Finally, it may be noted from Figure 9 that a phase shift 
exists between the position of the upper end of the cable and 
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the tension. Specifically, the tension tends to lag the 
displacement somewhat. This lag is caused by the dynamics 
of the phenomena; a juxtaposition of static solutions would 
necessarily show the tension and displacement to be in phase. 

The amount of computer time required to obtain numerical 
results depends on two factors, the number of subdivisions 
along the cable and the slope of the C g characteristic. In 
general, the number of calculations required, and consequent- 
ly, the computer time, is directly proportional to the square 
of the number of subdivisions along the s- axis or the cable 
length. The size of the time step advanced in each computa- 
tion is related to the slope of the C g characteristic 
according to Dt = Ds/C g where Ds is fixed by selection of 

the number of subdivisions of the cable and C is defined as 

s 

w 

C = (K/M) . Thus, in order to reduce the required computer 

b 

time it is necessary to use the minimum number of subdivisions 
of the cable that will give accurate results and use small 
values of C . Unfortunately C is dependent on the cable 
physical properties and, accordingly, is fixed. 

In view of the definition of C g it is apparent that the 
time steps for very "stiff" cables with large elastic con- 
stants K will be small. However, for lighter cables with 
smaller elastic constants such as would correspond to nylon 
rope, the time steps will be larger and the calculation 
procedure more efficient. For the particular example worked, 
30 subdivisions on the cable were used and it took approxi- 
mately 30 minutes of computer time to compute 2.9 cycles 
corresponding to 29 seconds in real time. 
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The value of C = (K/M) 2 has the physical significance 
s 

of being the speed of propagation of a stress wave along the 
cable while the second characteristic = (T'/M') 2 has 

the physical significance of being the speed of propagation 
of a transverse wave along the cable. Accordingly, for all 
physically realistic cables the stress wave speed would 
normally be much larger than the transverse wave speed and 
consequently, C g would be greater than C n » In the case of 
steel cables, the stress wave speed is normally quite large 
making the time step small. Moreover, in the limiting case 
of the inextensible cable the calculation procedure breaks 
down because the time step becomes zero. However, in this 
case or in cases where the cable stretching can be disregarded, 
the problem actually becomes much simpler. That is, the 
system of four equations (57-60) are replaced by the two 
equations (59) and (60) which are valid along the ±C n charac- 
teristics only. The ±C g characteristics are lost completely 
and the system reduces to a much simpler two-characteristic 
system. 
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V. CONCLUSIONS 



The mathematical model presented in this thesis provides 
a method for the solution of cable dynamic problems. Although 
a simplified boundary condition was employed on the surface 
for purposes of generating numerical results, the results 
obtained indicate that the method is valid and can be used in 
connection with any type of boundary condition on the surface. 

Excessive computer time is the greatest problem to be 
overcome in producing numerical results and two methods for 
reducing these were discussed. In cases where stretching can 
be disregarded the computer time can be reduced by orders of 
magnitude by utilizing the two-characteristic model since 
the time step is in this case defined by Dt = Ds/C n and 
is normally small compared to C g . In cases where the cable 
is relatively stiff and yet long enough such that stretching 
cannot be disregarded, no solution to the problem of excessive 
computer time presents itself. 
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RECOMMENDATIONS 



The following recommendations are advanced on the basis 

of experience gained during the course of the present study: 

1. A surface boundary condition should be applied that 
contains all the terms involved in a surface vessel's 
motion including mass, added mass, damping and buoyant 
force. 

2. The two-characteristic model should be developed and 
compared to the solutions using the four-characteristic 
method in order to evaluate the effect of stretching. 

3. A study of the effect of cable subdivision size on 
accuracy should be carried out. 

4. The model appearing in this thesis should be extended to 
include tangential drag, and surface wave drag on the 
cable and should be extended to the three dimensional 
case. 
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APPENDIX B 



It is interesting to note that if the current is zero 
then equation (63) becomes: 



dT _ sin l df 
T cos ¥ 



( 78 ) 



This equation can be easily integrated into the following form: 

cos H* 

T = (79) 

cos T 



By setting the current equal to zero in equation (64) it is 
not apparent that the result will be the same as equation (79) . 
It is possible to show that equation (64) does reduce to 
equation (79) by the following method. 



Letting 



£ = 



4EV 



and using a binomial expansion for (1+e) the following 
is obtained: 



(1+zft = 1 + je - i e 2 + yg- e 3 (80) 

Substituting equation (80) into the numerator of equation (64) 
produces the following numerator: 

, i 1 h u , 12 1 3 W 1 h „ . 1 12 1 3. 

(- 2 - ^ e cos e " T6 E ] ( " 2 G COS T + T ~ 8 E + 16 e ] 

(81) 
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The coefficient of the e 2 term will produce cos ¥ since the 
other terms will be of higher order. A similar analysis can 
be done in the denominator producing a coefficient of cos t . 
Equation (64) now reduces to: 



cos Y 

T = vT 

cos y 

which is identical to equation (79) . 



(82) 
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